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Abstract, fn this paper we present a method for conducting a coherent search 
for single spin compact binary coalescences in gravitational wave data and 
compare this search to the existing coincidence method for single spin searches. 
We propose a method to characterize the regions of the parameter space where the 
single spin search, both coincident and coherent, will increase detection efficiency 
over the existing non-precessing search. We also show example results of the 
coherent search on a stretch of data from LIGO's fourth science run but note that 
a set of signal based vetoes will be needed before this search can be run to try to 
make detections. 



1. Introduction 

The Laser Interferometer Gravitational-wave Observatory (LIGO) and Virgo scientific 
collaborations have performed many searches for compact binary coalescence (CBC) 
signals in data taken by gravitational wave interferometers pj [3] . The majority of 
these searches have utilized template waveforms where the spins of the individual 
components are neglected. In some areas of the parameter space spin can have 
a significant effect on the evolution of the system, and consequently the emitted 
gravitational waveform [U[5], leading to a poor match with the non-spinning templates. 
In these regions of parameter space the use of templates incorporating spin will provide 
an increase in search sensitivity. 

Incorporating spin into template waveforms in a gravitational wave search is a 
complex problem. In a non-spinning search for CBCs with circular orbits, a source 
is described by nine physical parameters [BJ. The majority of these do not affect the 
signal morphology but serve to change the overall amplitude, phase or coalescence 
time of the signal and are easily maximized over [7] . Therefore, template placement 
can be restricted to the two dimensional space of component masses [5] . A spinning 
CBC in a circular orbit, however, is described by 15 physical parameters [5J. The 
challenge is to formulate a method to detect any manner of spinning system while 
limiting the number of templates, such that an analysis can be run in a reasonable 
amount of time. The problem is simplified if the spins are aligned with the orbital 
angular momentum. In this case the system will have no precession and is described 
by just two extra parameters — the spins of both bodies in the direction of the 
angular momentum. Furthermore, these non-precessing waveforms are well described 
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by a single spin parameter [9], and it it therefore feasible to search for non-precessing 
waveforms using a three dimensional template bank. 

At the time of writing only one search for CBCs using spinning templates with 
precession in LIGO/ Virgo data has been published [TO]. This search utilized a 
phenomenological waveform family designed to capture precessional effects (TTJ, but 
was later abandoned because it was not found to increase efficiency relative to the non- 
spinning search [3j I12j . This was due to the ability of the phenomenological templates 
to match non-stationarities in the data and the lack of an effective signal consistency 
test to veto them such as the x 2 test used in the non spinning search [T5] . 

The physical template family (PTF) waveforms proposed in and further 
explored in [14l [15j [16] give a different method for searching for spinning binaries 
with precession. This method uses single-spin precessing waveforms as templates. 
Making clever use of maximization, it was shown [TJ] that a PTF search could be 
performed with a four dimensional template bank: the two masses, the magnitude 
of the spin and the angle between the spin and the orbital plane. This method is 
especially useful for detecting neutron star-black hole binary (NSBH) systems, where 
the spin of the neutron star would have a negligible effect on the dynamics of the 
system [TJ]. A coincidence search utilizing the PTF waveforms has been developed 
[16] . Data from each instrument is analysed separately and only events observed with 
consistent time of arrival, mass and spin parameters in more than one detector are 
retained. While coincidence requirements for non-spinning searches are well known 
[17] , it is less clear how to define coincidence when the additional spin parameters are 
present. 

In a coherent search [T51 [TU [501 |5TJ [55] the data from all active detectors are 
combined together before searching for interesting events in the combined data. This 
circumvents the need for a coincidence test between events in different detectors. 
Furthermore, a coherent search offers an increased detection efficiency over the 
coincident technique when more than two detectors are active |22| . The coherent 
technique is especially useful when the sky position is known, such as when searching 
for gravitational waves in coincidence with an electromagnetic transient, such as a 
gamma-ray burst (GRB) [53]. Since NSBH and binary neutron star (BNS) mergers 
are the preferred progenitor model [21] for short GRB, a coherent single-spin search 
is ideally suited to this source. 

In this paper we describe the implementation of a coherent search for single spin 
binaries with known sky location, using the PTF waveforms. We briefly review the 
PTF formalism before deriving the coherent PTF signal-to-noise ratio (SNR). Due to 
the increased complexity of the spinning waveforms, the coherent SNR has a different 
distribution than its non-spinning counterpart. In particular, there is a greater chance 
of obtaining a large value of the spinning SNR, even in Gaussian noise. Thus there 
is a trade-off between the improved spinning signal model and the increased false 
alarm rate at a fixed SNR. We explore the single-spin CBC parameter space to 
identify regions where spin (and precession) effects are significant enough to make 
the spinning search worthwhile. We will also briefly discuss some possibilities for 
vetoing background non-Gaussian transients in the data when using the PTF search 
and present results of this search run on a short stretch of LIGO's fourth science 
run (S4) data 

The layout of this paper is as follows: In section [2] we briefly review the 
single detector PTF search and investigate the distribution of the spinning SNR in 
Gaussian noise. In section [3] we introduce the coherent PTF search and investigate 
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the distribution of the coherent spinning SNR. In section [2] we identify regions of 
the parameter space where the PTF search offers increased sensitivity over the non- 
precessing search. Section [5] briefly describes our search pipeline and the results of 
these methods applied to a stretch of data from S4. 



2. Spinning Search Using Physical Template Family Waveforms 

In this section we give a brief recap of the PTF search and its implementation. We also 
explore the expected distribution of the spinning SNR in Gaussian noise and compare 
this to that of the non-precessing search. For a more detailed description of the PTF 
search and terminology we refer the reader to [TU [16] . We will follow the conventions 
of these earlier publications as much as possible. We will also assume that the reader 
is familiar with matched-filtering techniques and its application to gravitational wave 
data analysis, if not we refer the reader to [2"5l 17]. 



2.1. Single detector analysis 

The likelihood ratio of there being a signal h present in the data s for a single detector 
is given by 

\ / j \ P{s\h) 

Assuming the noise is Gaussian, the log likelihood can be written as 

logA=(/i| S )-i(/ l |/i). (2) 
Where we have defined the single detector inner product 

and Sh(f) is the power spectral density (PSD) of the detector. From this starting 
point, h must be re-expressed in such a way that it is possible to maximize over the 
majority of the parameters, leaving us with only a small number of dimensions over 
which to carry out a templated search. 

The dominant harmonic of the gravitational waveform can be expanded in terms 
of the five I = 2 spin-weighted spherical harmonics. The amplitude of each of these 
terms will depend upon the distance to the source, D; the sky location of the source, 
(9, ip); the orientation of the source, which is described by three angles: the inclination, 
l, polarization, <f>, and the orientation of the spin in the orbital plane, Lp. The waveform 
for each of these harmonics depends on the two masses, (M\,M.2); the amplitude of 
the spin, x~> t ne angle between the spin and the orbital plane, k; the inital orbital 
phase relative to the spin direction, $ 0j an d the time of coalescence, t c . Consequently 
the gravitational waveform for a single spin binary can be expressed as 

5 

h (t) = ]T Pi{D, 9, 0, l, l p)Q I (M u M 2 , X , «, $o, t c ). (4) 
i=i 

where Pj are five amplitudes and Q 1 describe five waveform components. 

To obtain the PTF detection statistic, a free maximization is carried out over the 
five Pi components as well as the initial orbital phase of the system. The SNR will then 
depend on 10 components: the and f phases of the five Q 1 waveforms. Specifically, 
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we calculate the inner product between each component Q 1 of the waveform and the 
data 

A 1 =(3^1) and S J = (s|Qi), (5) 
as well as between the different Q 1 themselve^J] 

M IJ = {Ql\Qi) = {Q\\Qi). (6) 
The maximized PTF detection statistic is given by [TH [T5] 



p , . i 1 

— = logA| max = - 



A T M _1 A + B T M *B 



+ i y / (A t M" 1 A - B T M-!B) 2 + (2A T M- 1 B) 2 .(7) 

The expression for the SNR, p, can be simplified by performing a transformation 
such that both Qq and Qi are orthonormal. First, perform a rotation on the Ql to 

make M IJ diagonal, then normalize the basis vectors. We denote the orthonormal 
basis Qq. This transformation will also orthonormalize Qi and render M IJ the 
identity matrix. After this transformation, the SNR can be written as 



p l = (^A-A + B-B) + MA- A- B ■ B) 2 + (zA-B^ . (8) 

where A and B are defined as in (j5j). 

We have performed a free maximization over the five Pi amplitudes. In principle, 
these depend upon six physical parameters. However, these parameters only enter in 
four different combinations as 

• an amplitude parameter, dependent on (D,9,^,cj)), 

• the relative sensitivity of the instrument to the + and x polarizations, dependent 
on (6, ip, 4>) 

• the inclination angle, t 

• the spin orientation, ip. 

Therefore performing a free maximization over the five Pj components means that 
the maximized Pj values may not correspond to a physical set of parameters. This is 
discussed in |14j and various methods for projecting onto the physical sub-space have 
been proposed. For the case of an externally triggered search, where the sky location 
is known, the situation is unchanged as the Pj are still described by the same four 
unknown parameters. 

When the orbital plane of the system does not precess, there is gravitational wave 
emission in only two of the harmonics, Q 1 and Q 2 . The other components vanish 
identically. Furthermore, these two harmonics are related by a phase shift: Q 1 = iQ 2 . 
Thus, the matrix M is degenerate and the PTF maximization breaks down. It is, 
however, straightforward to maximize over the two remaining amplitudes, and obtain 
the SNR as 

2 _ {s\Ql) 2 + {s\Q\) 2 

9 (QhlQh) i} 

t We have made the standard assumption that Qq = iQ\ . 
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This is identical to the well known SNR for the non-spinning search [7], and the two 
phases of Q 1 correspond to the and ir/2 phases of the non-precessing template. 

The PTF search allows one to perform a search using single spin waveforms in a 
reasonable amount of time on a single detector [IB]. Any event with an SNR above 
some preset threshold constitutes a single detector "trigger", and candidate events 
would be required to be observed in more than one detector. However, it remains 
a challenge to derive a metric on the four dimensional mass and spin space that 
could be used in generating a template bank and in defining coincidence requirements. 
Furthermore, a strategy for vetoing non-transient glitches has been suggested [TBJ, such 
a strategy would be needed to make a coincident PTF search viable. 

2.2. SNR distribution in Gaussian noise 

The PTF template waveform will provide a better match than a non-spinning template 
to a gravitational wave signal from a spinning binary. However, we pay a price since 
we must filter the data against more waveform components, Q 1 , thereby increasing 
the chance of a spurious match with the noise. Additionally, the spinning SNR takes 
a more complex form Q than the simple quadratic expression ^ when there is no 
precession. Here, we will investigate the SNR distributions in Gaussian noise for these 
two cases. In section|4]we use this to identify regions of parameter space with sufficient 
spin effects to warrant the use of the PTF search. 

Ten filters are used in the calculation of the PTF detection statistic: (Qq, s) and 
(Q K ,s). As both Qq and Qi are orthonormal, the only remaining freedom is the 
relation between the Qi, and Qi terms, 

N IJ = {Ql,Ql). (10) 

This N IJ is a will be a 5 x 5 antisymmetric matrix which can have 4 non-zero 
eigenvalues: ±Ai,±A2. The values of these eigenvalues determine the distribution 
of the PTF detection statistic. 

For every NSBH waveform we have tested using the initial LIGO sensitivity 
curve, the magnitudes of Ai and A2 have been very close to 1. Thus, although there 
are ten different waveform components, we find that, in effect, only six of these are 
independent — the others are linear combinations of these six. There are then only six 
independent filters and it is not difficult to show that the spinning SNR ^ collapses 
to a quadratic form which is x 2 distributed with six degrees of freedom in Gaussian 
noise. This is the "best" case for the detection statistic. The "worst" case occurs 
when both Ai and A2 are zero and all ten of the filters are independent. In this case, 
the SNR expression cannot be simplified and its distribution does not correspond to a 
X 2 distribution with 10 degrees of freedom as might be expected; the real distribution 
is somewhat more complex. Both best and worst cases are illustrated in Figure [l] 

The SNR (|9| for a non-precessing template follows a \ 2 distribution with two 
degrees of freedom in Gaussian noise. This is also plotted on Figure [T] By comparing 
the distributions of the PTF and non-precessing SNRs, it is clear that the background 
triggers produced by the PTF search will have, on the average, considerably larger 
SNR than those produced by the non-precessing search. We explore the effect that 
this has on a search further in section [4] 
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Figure 1. The "best" and "worst" possible distributions of the single detector 
PTF SNR squared, this is compared with the non precessing SNR squared. 



3. Coherent Spinning PTF search 

In this section, we introduce a multi-detector, coherent formulation of the PTF 
search. As in [55], we will restrict attention to the case where the sky location is 
known. This simplifies a coherent search as the sensitivity of the detectors to the 
two GW polarizations and the relative time delays between detectors are known. 
Astrophysically, this is of interest when searching for gravitational waves associated 
to electromagnetic transients such as GRBs [2H [26] . 

The data from various detectors are combined together coherently to form 
two coherent data streams, with one stream containing the + polarization of any 
gravitational wave signal present in the network and the second containing the x 
polarization. The coherent method will offer an improvement in sensitivity over the 
coincidence method when more than two detectors are used, as only two data streams 
are searched. For networks with greater than two detectors, it is also possible to 
construct null streams which will contain no gravitational wave signal, and can be 
used as a consistency test [371 H2] • 

We begin by formulating the coherent SNR for the spinning PTF search and go 
on to explore how this will be distributed in Gaussian noise. 



3.1. Coherent PTF Search Method 

To formulate a coherent detection statistic for the PTF templates we draw on many 
of the methods and techniques that were used in deriving the single detector statistic 



described in section 2.1 and derived in detail in [THITB]. We follow the conventions of 
[22] in extending this to a coherent, multi-detector search. Assuming that the noise 
in different detectors is independent, the multi-detector likelihood is given by 

kA = (h|s)-~(h|h) (11) 
where we have defined a multiple detector matched filter 



(a|b) = 5>*|6*). (12) 

x 

and the index X runs over the detectors in the network. 
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As before, we want to maximize over as many of the parameters as possible to 
minimize the dimension of the required template bank. We start by maximizing this 
over the distance, D, and initial orbital phase, $o, to obtairj|] 

\2 , / „ y t-, y \ 2 



lnA| 



max(D,<E> ) 



1 Ex,/ 
2 



(13) 



where Pj are the amplitudes of the various waveform components Q 1 , and A^, B J X 
are defined as in ([5]). Although the Pj depend upon D the maximized likelihood is 
independent of it as scaling the distance has an identical effect on both the numerator 
and denominator of eq. (13). 

As in the single detector case, we would like to maximize over the Pj to eliminate 
them. However, in the mult i- detect or case, they are detector dependent since the 
sensitivity of the detectors to the + and x gravitational wave polarizations will differ. 
These sensitivities are encoded in the detector response functions, F + and F x , which 
depend on the sky location of the source in the detector frame. As we are focusing on 
an externally triggered search, where the sky location is known, these values will be 
known for each detector. We can then factor the detector dependent terms out of the 
P? as 

Pf = F*{0, 4>)Si{D, l, t(,, tpo) + F*(fi, <f>)T T (D, l, J>, <po) (14) 

where Si and 7/ denote the amplitude of the + and x components respectively of 
the 5 Qi in the radiation frame. They depend on the distance, D and the angles 
(L,<p,tpo) that describe the rotations necessary to transform from the source frame to 
the radiation frame. 

We can re-cast the log-likelihood into a form which more closely resembles the 
single detector case by introducing ten-dimensional analogues of the Pi and Q 1 by 
defining 

V a '■= [Si, S%, S3, S4, S5, Ti, T2,Ts, T4,, Tg] 



Go. 



(15) 



: F +Qo,% ; ■■■;^+ ( 3o,f ;F X Q * ; ...;P X <9 

The change to ten dimensions naturally arises because a multiple detector coherent 
network is sensitive to both the + and x components, whereas a single detector 
network is only sensitive to one polarization. We also define the multi-detector inner 
products between signal and waveform components 

A a =(810?) and B a = (s\Q%) 



M ap = (esieg) = (efie|. 

The log likelihood equation can then be written as 

1 V a Vp (A a A 13 + B a BP) 
2~ 



InA], 



<D,* ) 



(16) 



We proceed, as before, by transforming to an orthonormal basis QqiG" f° r the 
waveform components. Then, maximizing freely over V a yields the coherent PTF 
SNR 



2 

Pcoh 



A-A + B-B 



A-A-B-B 



2A-B 



(17) 



§ The Y subscript in the inner product in the denominator denotes the fact that the PSD of detector 
Y is used in evaluating the inner product. We do not require the noise PSDs of the different detectors 
to be the same 
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Figure 2. The "best" and "worst" possible distributions of the coherent PTF 
SNR squared as well as the distribution of the non spinning SNR squared. 



where, as before, the tilde denotes that we are in the orthonormal basis. This is very 
similar in form to the single detector statistic in equation Q. When the network 
is only sensitive to one polarization, the matrix A4 a " becomes degenerate and the 
maximization procedure must be re-visited. Here it is natural to remove all terms 
corresponding to the second polarization and reduce to 5 dimensions, as in the single 
detection search. Additionally, in section [2~l] we noted that when the template has no 
precession the single detection PTF SNR collapses to the familiar SNR formalism used 
in the non-spinning search. Similarly, in the coherent PTF search, when the template 
has no precession, the coherent SNR will collapse to the non-spinning coherent SNR 
given in [2"2"| . 

The coherent SNR of equation (17) can be used as a detection statistic in 
performing a coherent search using PTF templates, as we explore in section [5j In the 
single detector search, we maximized freely over five Pr which were dependent upon 
four physical parameters. Here, the V a still depend on only 4 parameters but we are 
now maximizing over ten amplitudes. This clearly introduces a lot of unnecessary 
freedom. We are currently investigating alternative methods of constructing the 
coherent SNR which might eliminate these un-physical degrees of freedom. However, 
we should note that the coincidence search allows for a similar freedom as the Pj are 
maximized independently for each detector. Consequently, for a network with three 
or more detectors, the coherent search provides a sensitivity improvement. 



3.2. SNR distribution in Gaussian noise 

In section |2.2| we explored how the single detector PTF statistic is distributed in 
Gaussian noise. For the coherent PTF search we can use a similar strategy to 
investigate the distribution of the coherent SNR. In the coherent case there are twenty 
filters A a and B a and we have constructed the detection statistic such that Qq and 
Q" are orthonormal. As before, the only freedom is the relationship between the 

2 

and ^ terms encoded in 

ti a0 = (Q o ,Qi). (18) 

This is a 10x10 antisymmetric matrix comprised of four 5x5 blocks, each of which 
is antisymmetric. Therefore this matrix can have 8 non-zero eigenvalues: iAi.,2,3,4. 
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Non Precessing 


PTF 


Single detector 


6.79 


7.63 


Coherent 


7.26 


8.53 



Table 1. The SNR corresponding to a FAP of lO -10 for the non precessing and 
the PTF search, for both coherent and single detector cases. Here, for the PTF 
case the single detector and coherent detection statistics are assumed to be \ 2 
distributions with 6 and 12 degrees of freedom respectively. 



These eigenvalues determine the distribution of coherent SNR in Gaussian noise — for 
smaller eigenvalues, the large SNR tail of the distribution becomes more significant. 
In the tests that we have performed using the initial LIGO sensitivity curve and NSBH 
precessing templates, all four eigenvalues give values close to unity, the "best" case 
in which there are 12 independent waveform components. However, the distribution 
does not collapse to a x 2 - In Figure [2] we demonstrate that this gives a distribution 
similar to a \ 2 distribution with 12 degrees of freedom. In the "worst" case, where 
all of the eigenvalues are equal to 0, there are 20 independent waveform components 
and this distribution is also shown in Figure [2] 

4. Identifying where the PTF search is most beneficial 

In sections [2] and [3] we have derived the spinning SNR that can be used to perform a 
gravitational wave search using single spin inspiral waveforms as templates. We have 
demonstrated that, on the average, background triggers will have larger values of SNR 
in the PTF search than in the non-precessing search. At the same time, precessing 
PTF waveforms will be a better match to any spinning, precessing signals in the 
data. This begs the question as to whether it is preferable to use a search with non- 
precessing waveforms or single spin PTF waveforms to detect precessing systems. The 
PTF triggers will match the waveform better but this comes at the cost of searching 
a larger parameter space. 

To quantify this, in Table [l] we give the SNR that corresponds to a false alarm 
probability (FAP) of 10~ 10 in Guassian noise for the various searches. We chose this 
value because it roughly corresponds to the loudest background events we observe 
when running the search on 2000 seconds of Gaussian noise, as is appropriate for a 
GRB search. The figures in the table show that the PTF search must obtain 26% 
more signal power (SNR squared) to be more efficient in the single detector case at 
this FAP and 38% more signal power for the coherent case. 

There are large areas of the parameter space where precession will not significantly 
effect the evolution of the binary and thus a non-precessing template will pick up the 
majority of the power in a precessing signal. In these areas it would be better to search 
for the spinning signal with a non-precessing template, achieving a lower FAP than for 
the PTF search using an exactly matching template. Equivalently, when a system has 
little precession, the majority of power is contained in the Q 1 and Q 2 components of 
the PTF waveform and these two components are very similar, up to an overall phase 
shift. We can then consider performing a "restricted PTF" search, where we filter 
only these two components of the waveform against the data. This serves to reduce 
the FAP at a fixed SNR while losing only a small amount of the power in the signal. 

To do this, we test every template waveform, before filtering, to determine 
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Figure 3. The fraction of templates analysed by the full PTF statistic as a 
function of the masses in the NSBH region of the parameter space. 

whether on whether the template would be more likely to detect a matching signal 
below a false alarm probability of 10" 10 using the restricted or full PTF search. This 
can be calculated by simulating a large number of gravitational wave signals, with 
masses and spin matching those of the template, but uniformly distributed in volume 
and orientation. Then, simply count the number of simulated signals expected to give 
an SNR greater than the value corresponding to a FAP 10 -10 (given in Table [TJ for 
both methods. Whichever of the PTF or restricted methods is expected to perform 
better is then used when filtering the data with that template. Using this method, 
we are able search the full parameter space of NSBH binaries in a single search, 
including non-spinning, non-precessing, marginally precessing and fully precessing 
configurations. This method works equally well for the single detector or the coherent 
search. 

In Figure [3] we illustrate the fraction of templates analysed by the full PTF 
statistic, as a function of the masses, for the coherent search. The splitting of the 
templates into full and restricted does not require filtering against the data, but it 
does make use of the PSDs of the detectors. For this study, we use data from the 
three LIGO detectors during the S4 run. 

A template bank was generating by taking a standard non-spinning template bank 
[28] in the mass space and, for each value of the masses, creating 15 templates with 
identical masses but spin parameters gridded over the two dimensional spin space, 
as described in [16] . The precessing single spin templates are most needed in the 
high mass ratio region of the parameter space. For this template bank, there are 
35395 templates to be analysed with the restricted method and 14660 templates to be 
analysed with the full PTF method. 

5. Search method and example results 

In section[3]we derived a detection statistic appropriate for a coherent search using the 
PTF waveforms as templates. In section [4] we described a method through which one 
can identify where the PTF search is most needed and to split a template bank into 
those templates that should be analysed with the full PTF statistic and those that 
should be analysed with the restricted PTF statistic. We have combined these two 
methods together to create a search pipeline than can be used to coherently analyse 
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gravitational wave data to search for precessing NSBH signals associated to short 
GRBs. We will briefly describe the analysis procedure before presenting an example 
result. 

The search uses much of the same architecture as that described in [22] and [2S] . 
Namely we search for gravitational wave signals in the "on-source" time, defined to be 
[-5,+l) seconds around the reported time of the GRB. Background is estimated from 
performing 324, 6 second trials around the GRB time, but at least 48s away from the 
on-source time. The coherent PTF search makes use of the same infrastructure as a 
coherent non-spinning search described in |22j . In particular, the data handling, PSD 
estimation and matched filtering routines are the same. Of course, the coherent PTF 
search makes use of spinning, precessing waveforms in the filtering and computes the 
SNR given in equation (17). 

To demonstrate the performance of the coherent PTF, we ran it over a stretch 
of data from LIGOs S4 run. The data was chosen randomly, subject to the condition 
that all three of the LIGO detectors were operational at the time. This is the same 
data as was used to illustrate the template bank splitting in section [4j and the same 
bank with 15,000 full PTF and 35,000 restricted PTF templates was used. 

In Figure [4] we show the distribution of the SNR of the triggers produced using 
both the full and restricted statistic. This is shown for the stretch of real S4 data 
and for a stretch of simulated Gaussian data. As expected, the SNRs of triggers in 
Gaussian noise are larger for the precessing templates than the restricted ones, even 
though significantly fewer templates were analysed with the full statistic. The results 
from real data are badly affected by non-Gaussianities in the data. A number of loud 
transients are clearly visible as short duration peaks of large SNR, while there are 
an even greater number of quieter peaks throughout the analyzed time. This has a 
similar effect on both the full and restricted waveforms. 

In |22j . we described and developed a number of tools which can be used to 
effectively remove the majority of the non-Gaussian features from a non-spinning, 
coherent analysis. These include null stream consistency [27 , amplitude consistency 
and x 2 signal consistency tests [131 [3D]. All of these can be applied without 
modification to the restricted PTF search, and it seems reasonable to expect they 
would be similarly effective in reducing the effect of non-Gaussianities in the data. 
However, we currently have no such tools which can be used for the full PTF 
waveforms. Before using this search on real data we will need to implement a set 
of tests that can discriminate glitches from real signals for the full statistic. It 
should be relatively straightforward to implement the null stream consistency test. 
Unfortunately, as discussed in [22 , the null stream for the LIGO S4 detectors is 
constructed only from the two instruments in Hanford. In this stretch of data the 
loudest background triggers are caused by non-stationarities in the Livingston detector 
and thus the null stream is ineffective. Alternatively a x 2 test such as the ones 
described in [T3J [3D] could be adapted to this search, [TB] presents a possible way of 
doing this for single detectors. We are working on developing an alternative version of 
this x 2 test, which would test the consistency of the six independent components of a 
single detector PTF waveform, and then extending this to the fully coherent analysis. 



6. Discussion 



In this paper we have presented a method for performing a coherent search for 
precessing, single spin black hole-neutron star coalescences using the PTF method. 
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Figure 4. The distribution of triggers found by the coherent PTF search. The 
left panels show the distribution of triggers from templates that were analysed 
using the "restricted" coherent PTF search, the right panels show the distribution 
from the templates that were analysed using the "full" coherent PTF search. The 
top panels were created from analysing Gaussian noise. The bottom panels were 
created from analysing a stretch of real data from S4. All these plots have been 
rescaled to use the same y-axis. For the two cases using real data the non-Gaussian 
spikes extend much higher than is shown, the loudest trigger has an SNR of 39 in 
the restricted case and 45 in the full case. 



We have compared the performance to searches using non-precessing waveforms and 
have identified regions of the parameter space where the PTF search offers increased 
sensitivity. We have presented a method by which these areas could be identified and 
demonstrated these techniques on a short stretch of S4 data. 

This method should allow for the detection of highly precessing NSBH sytems 
with greater efficiency than the current non-spinning searches. However, more work is 
required before this search is ready to be used. The main need is for the development 
of effective methods of separating glitches from real events in the full PTF search, 
whether performing a coincident or a coherent search. As discussed in section [5j it 
should be possible to adapt a lot of the methods that have proven effective in non- 
spinning searches [22] but this is a non-trivial task. 

In this paper, we have focused on the PTF precessing waveforms. However, many 
of the techniques we have discussed would be equally applicable to the dominant 
harmonic of any family of precessing waveforms. In particular, the method of 
maximizing over freely over the amplitudes of the five components of the I — 2 spin 
weighted-spherical harmonic is directly applicable to other waveform families. As the 
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catalogue of numerical simulations of precessing binaries grows, these methods may 
well find applications in searches using such as numerical relativity inspired inspiral- 
mcrger-ringdown waveforms. 
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